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A simple approximate expression in real and reciprocal spaces is given for the static exchange- 
correlation kernel of a uniform electron gas interacting with the long-range part only of the Coulomb 
interaction. This expression interpolates between the exact asymptotic behaviors of this kernel at 
small and large wave vectors which in turn requires, among other thing, information from the 
momentum distribution of the uniform electron gas with the same interaction that have been cal- 
culated in the Go Wo approximation. This exchange-correlation kernel as well as its complement 
analogue associated to the short-range part of the Coulomb interaction are more local than the 
Coulombic exchange-correlation kernel and constitute potential ingredients in approximations for 
recent adiabatic connection fluctuation-dissipation and/or density functional theory approaches of 
the electronic correlation problem based on a separate treatment of long-range and short-range 
interaction effects. 

PACS numbers: 71.10.Ca,71.45.Gm 



I. INTRODUCTION 

Density functional theory (DFT) 1 applied within the 
Kohn-Sham (KS)£ scheme is nowadays a widely-used 
method for electronic calculations in condensed-matter 
physics and quantum chemistry. The design of better 
approximations to the exchange-correlation energy func- 
tional, the central quantity of the theory, constitutes an 
important topic of research. 

Following early ideas^* 5 -, the possibility of improv- 
ing the description of the exchange and/or correlation 
energy functionals by treating separately its long-range 
and short-range components has recently gained grow- 
ing interest 6 . 7 ! 8 ! 9 ! 10 . 11 ! 12 . 13 ! 14 : 15 ! 16 ! 17 ! 18 ! 19 ! 20 ! 21 : 22 . In this 
approach, the Coulomb electron-electron interaction is 
decomposed as 



v»(r) +u^(r), 



(1) 



where v^ir) is a long-range interaction, v^ir) is the 
complement short-range interaction and /i is a pa- 
rameter controlling the range of the separation. For 
example, a convenient long-range interaction which 

has been often used is the so-called erf interac- 
ti3n B.fl.in.n.i2. 1 4.is. 1 fi. 1 7.i8. 1 fl.2n.2i.2a 



<rf( r ) 



erf(/xr) 



(2) 



which vanishes for fi = and reduces to the Coulomb 
interaction for /i — ► oo. The Coulombic exchange- 
correlation energy functional E xccou i[n) can be in turn 
decomposed aa 6 . 7 . 8 . 9 . 10 . 12 . 20 . 21 



in n \ 



EL 



,erfl 



(3) 



where E^ c erf [n] is the long-range component of exchange- 
correlation energy, associated to the interaction w^ rf (r), 
and E£ c eli [n] is the complement short-range part. 



The long-range exchange-correlation energy, E£ crf , 
can be efficiently calculated by configuration interac- 
tion (CI)2ii£ or multi-configurational self-consistent-field 
(MCSCF) 2 ^ methods enabling to handle near degener- 
acy in atoms and molecules, and by second-order per- 
turbation theory 25 or adiabatic connection fluctuation- 
dissipation (ACFD) 13,26 approaches for describing van 
der Waals interactions. The short-range exchange- 
correlation energy E*t f can be well approximated by 
simple (semi) local density functional approximations^!^ 
based a uniform electron gas model with a modified 
electron-electron interaction&S 7 . (see also Ref . l28j) . 

Knowledge of the properties of the uniform electron gas 
is always useful to gain more insight into the physics of 
electronic correlation of inhomogeneous systems and ulti- 
mately to construct better approximations. In particular, 
in standard KS-DFT, the static exchange-correlation ker- 
nel of the Coulombic uniform electron gas, f X c,cou\, is of 
interest and several accurate parametrizations are avail- 
able (see, e.g., Refs. I2gl3fll3l|) . Indeed, f xc , C oui has been 
recognized as a potential ingredient for density functional 
approximations (see, e.g., Refs. Hl l2l32l3 3 34 35 3a l37|) . 
fxc,coui can also be used in principle in ACFD approaches 
(see, e.g., Ref.l38lfor an example using a simplified kernel) 
or time-dependent density functional theory (TDDFT) 
for density-density response calculations, although in 
these cases clear strategies to use f X c,cou\ are still under 
investigations. As a matter of fact, the high non-locality 
of fxc,coui questions its transferability to inhomogeneous 
systems and limits its use. 

In the framework of a long-range/short-range energy 
decomposition in DFT, knowledge of the static exchange- 
correlation kernel fj* f of a uniform electron gas with 
the long-range erf interaction tr rf (r), as well as its com- 
plement analogue f£ Cjecf = /xc,coui - fxc,mxi associated 



to the short-range part of the Coulomb interaction, is 
also of interest. In particular, f£ crf can be used in ap- 
proximations for the long-range exchange-correlation en- 



2 



ergy E^ c erf within ACDF approaches and f£ c erf can be 
used to construct density functional approximations for 

fol- 



,orf ' 



the short-range exchange-correlation energy E£ c 
lowing the same ideas as for the Coulombic functional 
E xc ,cou\- In comparison to / X c,coul> the non-locality char- 
acter of f£ c erf and f£ c erf tends to be attenuated by the 
reduction of the interaction and the modified kernels may 
thus be more transferable to inhomogeneous systems. In 
principle, these kernels can also be useful in the context 
of TDDFT. 

The main purpose of this work is therefore the deter- 
mination of the static exchange-correlation kernel of the 
uniform electron gas with the long-range erf interaction 
/xc,erf> the complement kernel /£ c erf = / XC)COU i - fg c>(ai 
being directly deducible from it. A simple approxima- 
tion is constructed for fj* crf from its exact asymptotic 
behaviors which in turn require, among other things, in- 
formation from the momentum distribution associated to 
the modified interaction. Therefore, the momentum dis- 
tribution is preliminary calculated by many-body per- 
turbation theory (MBPT), in particular using the GW 
approximation**. All the calculations are also performed 
with the Coulomb interaction which enables a check by 
comparison to other available data. The production of 
data of the same quality for /a; Cl coui and f£ c orf is also 



rf 

important to construct the complement kernel f xc erf 



Let's describe now the system under consideration. 
The uniform electron gas is described as N electrons in a 
box of volume Q with a uniform neutralizing background, 
studied in the thermodynamic limit (i.e. N — > oo and 
£1 — > oo such that the density n = N/ft remains con- 
stant). In second quantization, its Hamiltonian writes 

H = f + V, (4) 
where T is the kinetic energy operator 

T = Y,< k )<i° c ^ ( 5 ) 

ka 

with e(k) = k 2 /2, and V is the electron-electron interac- 
tion 

^ = 2fi /ZYl U (9) C k+q<T C L-q<T' C k'^C k(T , (6) 

where the constant term q = has been omitted 
since it cancels out with the electron-background and 
background-background interactions. In Eq. ©, v(q) is 
the Fourier transform of an arbitrary electron-electron 
interaction. The standard case corresponds to the 
Coulomb interaction 



47T 

Wcoui(g) = -j, 



(7) 



while for the long-range erf interaction the Fourier trans- 
form writes 



V erf(>l) 



The paper is organized as follows. In Sec. [H] we dis- 
cuss the calculation of the momentum distribution and 
associated properties. In Sec. IIIII the static exchange- 
correlation kernel is obtained by interpolation from its 
limiting behaviors. In Sec. IIVI as a basic example of the 
use of this kernel, the long-range correlation energy of the 
uniform electron gas is calculated. Sec. summarizes 
and concludes this work. Details of derivations and an- 
alytical parametrizations can be found in Appendices lAl 

to El 

Unless otherwise stated, atomic units is assumed 
throughout this work. 



II. MOMENTUM DISTRIBUTION 

One of the most reliable calculations of the momen- 
tum distribution of the Coulombic electron gas for a wide 
range of densities are that of Takada and Yasuhara 40 
using the effective-potential expansion method, and the 
quantum Monte Carlo (QMC) calculation of Ortiz and 
Balloneii. Here, we will consider the use of the more 
traditional methods of MBPT. 

In MBPT, the momentum distribution n(k) is ex- 
pressed as 



i(k) = 



+ 00 



^ 0+ G(fc,.) 



(9) 



where G(k,uj) is the one-particle Green function, calcu- 
lated in a given approximation. Daniel and Voskoa* have 
calculated the momentum distribution of the Coulombic 
electron gas in the RPA. The fact that the RPA prob- 
ably constitutes the simplest MBPT approximation for 
the electron gas and becomes exact in the high-density 
limit (r s — ► where r s — (47m/3) -1 / 3 ) makes it valu- 
able for comparison to more elaborate approximations. 
Later, additional exchange terms have been included by 
Lam 4 ^, but the improvement over the RPA remains mod- 
est. More recently, more accurate MBPT calculations 
of the momentum distribution of the Coulombic electron 
gas have been reported using the GW approximation^ at 
different self-consistency levels^ 4 ^*^*^. In its semi-self- 
consistent, GWq, or fully self-consistent, GW, versions 
this method is known to be particle-conserving 45-48 - 49 , 
i.e. fulfilling the exact condition 



i(k) = 



(10) 



while the non-self-consistent GqWq approximation is 
However, even in this case, the violation of 



no 



-50 



(8) 



Eq. HlOfl is numerically small (2% at most for r s < 10 ac- 
cording to our own calculations). Moreover, imposition of 
the self-consistency does not necessary improve the mo- 
mentum distribution; in particular the self-consistency 
increases the quasiparticle renormalization factor at the 
Fermi surface Zp [Eq. Ill 111 ^ which is suspected to be 
already too large at the GqWo level (see below). We 



3 



therefore choose to compute the momentum distribution 
of the uniform electron gas with the modified erf inter- 
action in the Go Wo approximation (see Refs. I44l4fil5ll 
for details). For comparison purpose, calculati ons of the 
momentum distribution in the RPA (see Refs. I42l5ll for 
details) have also been performed. The range of densi- 
ties explored goes from r s =0.1 to r s = 10 and that of 
interaction parameters from fi = 0.1 to fx = 30. 

In order to assess the quality of the obtained momen- 
tum distributions, we look especially at two quantities 
that characterizes well the effect of correlations on the 
momentum distribution: the quasiparticle renormaliza- 
tion factor at the Fermi surface 



r 3 = 5. 



n(k 



(11) 



where kp — (ii^n) 1 ^ is the Fermi wave-vector, and the 
fraction of electrons in the correlation tail (i.e., k > kp) 



AN 



dk 
(2^ 



■Airk 2 n{k). 



(12) 



We also compute some moments of the momentum dis- 
tribution: the second reduced moment 82 related to the 
kinetic energy, and the fourth reduced moment $4. The 
dimensionless quantity 82 is defined as 



80 = 



(E K ) ~ {E K ) 



(E 



(13) 



K 



where Ep — (1/A) Sko- £ (^) c L (r Ck <' ^ s ^ ne kinetic energy 
per particle operator, (•••)o arid (•••) mean expecta- 
tion values between the non-interacting and interacting 
ground states, respectively. Similarly, $4 is defined as 



84 = 



(E%) (E 2 K )o 
(El), ' 



(14) 



where E\ = (1/A) £] ko . £(fc) 2 Ck CT c kcr . It is these mo- 
ments that appear in the asymptotic expansion of the 
static exchange-correlation kernel (see Sec. IHIj) . 

We first present results for the Coulombic uniform elec- 
tron gas and compare to available results in the literature 
in order to assess the accuracy of the Go Wo approxima- 
tion. 

Fig.nreports the momentum distribution n cou i(fc) cal- 
culated at the RPA and G W levels for r s = 5. The 
parametrization of Gori-Giorgi and Ziesche (GZ)£2 using 
accurate data from Takada and Yasuhara— and known 
theoretical constraints is also reported. One sees that 
the GqWq calculation greatly improves the RPA and is 
very closed to the GZ parametrization, especially for the 
correlation tail. 

The renormalization factor ^F.coui is plotted in Fig. [21 
Zf,com\ given by the RPA is largely underestimated and 
even negative at very low densities (r s > 7). As 
announced before, GqWq gives a renormalization fac- 
tor slightly too large compared to that of the GZ 
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FIG. 1: Momentum distribution with the Coulomb interac- 
tion for r s = 5: GZ parametrization— (solid curve), RPA 
(short-dashed curve) and Go Wo (long-dashed curve) calcula- 
tions. 
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FIG. 2: Renormalization factor with the Coulomb interac- 
tion: GZ parametrization— (solid curve), RPA (short-dashed 
curve) and Go Wo (long-dashed curve) calculations. 



parametrization which uses supposedly more accurate 
data of Takada and Yasuhara— for Zp :Con \. 

The fraction of electrons in the correlation tail 
(AAT/./V) C oui is reported in Fig. The RPA largely 
overestimates (AN/N) cou \. The G0W0 approximation 
gives (AA/A) cou i that is overall close to that given by 
the GZ parametrization, slightly deviating for large r s 
(r s > 8). As an independent check, we also report in 
Fig. |21 the fraction of electrons in the correlation tail ex- 
tracted from a Coupled-Cluster calculation with double 
excitations (CCD)— In this calculation, the energy of a 
uniform electron with unoccupied non-interacting levels 
shifted by a common gap is calculated and (AN/N) cou \ 
is deduced as the derivative of the energy at zero gap 
according to the Hellmann-Feynman theorem. The CCD 
calculation is close to both the GZ and Go Wo results. 

We now discuss results for the moments of the momen- 
tum distribution. 

Using the virial theorem, the reduced moment <52,coui 
for the Coulomb interaction can be expressed exactly 
in term of the correlation energy of the uniform elec- 
tron gas (see Appendix ^J. This enables us to per- 
form another test of the accuracy of the Go Wo calcu- 
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FIG. 3: Fraction of electrons in the correlation tail of the 
moment distribution with the Coulomb interaction: GZ 
parametrization^- (solid curve), RPA (short-dashed curve), 
Go Wo (long-dashed curve) and CCD^ (dotted curve) calcu- 
lations. 
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FIG. 5: (54. C o U i with the Coulomb interaction: GZ— (solid 
curve) and FHERp- (dotted curve) parametrizations, RPA 
(short-dashed curve) and Go Wo (long-dashed curve) calcula- 
tions. 
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FIG. 4: <52, C oui with the Coulomb interaction: formula from 
virial theorem [Eq. 1A5L thick solid curve], GZ parametriza- 
tion— (solid curve superimposed with the thick solid curve), 
RPA (short-dashed curve) and Go Wo (long-dashed curve) cal- 
culations. 
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FIG. 6: Momentum distribution with the erf interaction for 
r s — 5 and fi = 0.3: RPA (short-dashed curve) and Go Wo 
(long-dashed curve) calculations. For comparison, the Go Wo 
calculation for the Coulomb case (solid curve) is also shown. 



lation. Fig. 0] compares t>2 lC oui obtained with the virial 
theorem [Eq. (|A5I) ] using the usual Vosko-Wilk-Nusair 
(VWN) parametrization of the correlation energy^, with 
the GZ parametrization, and with the RPA and Go Wo 
calculations. The RPA calculation gives very poor re- 
sults for this quantity, especially for large r s . The GZ 
parametrization has been designed to reproduce the cor- 
rect kinetic energy and consequently gives <52,coui identi- 
cal to the formula from the virial theorem. The Go Wo 
approximation agrees well with this result for all r s . 

Fig. [5] shows the reduced moment J^coui given by the 
GZ parametrization and by the RPA and Go Wo calcula- 
tions. For comparison, a parametrization of <54 jCOU i due 
to Farid, Heine, Engel and Robertson (FHER)S is also 
reported. The G0W0 calculation agrees with the GZ 
parametrization for r s < 5 but deviates from above for 
larger values of r s . The FHER parametrization is con- 
stantly below the GZ result. These results reflect the 
fact that (J^coui is difficult to compute accurately due 
to its important sensitivity to the correlation tail of the 
momentum distribution. In comparison to the RPA cal- 
culation, the Gq Wq approximation, as well as the GZ and 



FHER parametrizations, give "reasonable estimates" of 
<54,coui- For later use, &4. cou i given by the GZ parametriza- 
tion and by the RPA and Go Wo calculations are fitted to 
analytical parametrizations in Appendix FBI 

Having assessed the precision of the Go Wo approxima- 
tion on the Coulombic electron gas, we now present the 
results for the long-range erf interaction. 

Fig. reports the momentum distribution n^ rf in the 
RPA and Go Wo approximations for r s = 5 and \x = 0.3. 
For comparison with Fig. the Go Wo calculation for 
the Coulomb interaction is also shown. The reduction of 
the interaction naturally brings the momentum distribu- 
tion closer to the non-interacting momentum distribution 
no(fc) = 9{kp — fc), even if the modifications are small at 
the scale of the plot, especially for the correlation tail. 
Also, the difference between RPA and Go Wo is reduced 
compared to the Coulombic case at the same density. 

The variation of the renormalization factor Zp f with 
the interaction parameter /i is reported in Fig. Q for the 
RPA and Go Wo approximation with r s — 5. One sees 
that Zp cif is very sensitive with respect to [i near /i = 0. 
Thus, even a very small interaction introduces significant 
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FIG. 7: Renormalization factor with the erf interaction for 
r s = 5: RPA (short-dashed curve) and GoWo (long-dashed 
curve) calculations. 



FIG. 9: &2 erf w ith the erf interaction for r s = 5: formula 
from virial theorem [thick solid curve, Eq. JXeJ], RPA (short- 
dashed curve) and Go Wo (long-dashed curve) calculations. 
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FIG. 8: Fraction of electrons in the correlation tail of the 
momentum distribution with the erf interaction: RPA (short- 
dashed curve), Go Wo (long-dashed curve) and CCD (dotted 
curve) calculations. 




FIG. 10: 5£ erf with the erf interaction for r s = 5: RPA (short- 
dashed curve) and Go Wo (long-dashed curve) calculations. 



incoherence into the evolution of an electron in the sys- 
tem. 

The fraction of electrons in the correlation tail 
(AAT/AOcrf is Plotted in Fig. for the RPA and G W P 
calculations. A CCD calculation similar of that of Ref.l53l 
but with the erf interaction is also reported. As for the 
Coulomb case, the Go Wo and CCD agree well in all the 
range of fi. 

The reduced moment 5% f associated to the long-range 
erf interaction can also be expressed in term of the corre- 
lation energy via a generalization of the virial theorem to 
this modified interaction (see Appendix |A"|) . Fig. com- 
pares <?2 er f derived from the virial theorem 
usin g t he parametrization of the correlation energy of 
Ref. |23 with the RPA and G W results. The G W cal- 
culation is in overall agreement with the formula from 
the virial theorem. The slight inaccuracy of Go Wo in 
the Coulombic limit (/x — > oo) is inherent to the method 
which neglects part of the short-range interactions. On 
the contrary, the inaccuracy for small /i is of numerical 
origins: when the interaction is too small, the magnitude 
of the correlation tail is lower than the precision of the 



method (which does not exceed 10 for the momentum 
distribution in the present implementation). 

Finally, the reduced moment S^ erl calculated in the 
RPA and Go Wo approximations are reported in Fig. 
for r s = 5. Analytical parametrizations of these results 
are given in Appendix^] The parametrization of 6% f in 
the Go Wo approximation will be used in the next section 
for the determination of the exchange-correlation kernel. 



III. STATIC EXCHANGE-CORRELATION 
KERNEL 

The dynamical exchange-correlation kernel f xc (k,uj), 
or equivalently the dynamical local-field factor G(/c, uj) = 
—v(k)f xc (k,u}), is related to the linear dynamical (in- 
teracting) response function x(k,oS) via the Dyson-like 
screening equation 

X(k, w) = XoO, v) + Xo(k, lo) [v(k) + f xc (k, u)] x(k, w), 

(15) 

where xo(fc,cj) is the non-interacting dynamical re- 
sponse function. The static exchange-correlation kernel 



6 



is fxc(k) = fxc(k,uj = 0). 

The study of the local-field 



factor of the Coulom- 



bic electron g£ 
(see, e.g., Refs . 


is has generated 


an 


abundant 


iterature 


29 30 55 56 ^j^^^m^^mA 


l67ll6M6Qll7fil7Th 


. The QMC simulations 30 69 give probably 



one of the most reliable static exchange-correlation kernel 
fxc.coui(k). A remarkable feature of the QMC results is 
that fxc.coui(k) does not have much structure, and can be 
essentially reproduced by combining its small-fc constant 
value for k < 2kp and its large-fc two-term asymptotic 
expansion for k > 2kp2&. For the case of the long-range 
erf interaction, accurate data are not available to check 
this property. However, it seems reasonable to assume 
that no additional structure appears in the exchange- 
correlation kernel when the interaction is reduced. As a 
first approximation, we will thus construct fj~ f (k) from 
its limiting behaviors for the case of the erf interaction 
too. 

The limiting behaviors of f xc {k) for both the Coulomb 
and long-range erf interactions can be cast into the same 



form. For k - 
has the limit 



0, the static exchange-correlation kernel 



4-7T 

fxc(0) = ~T2~A, 
Kp 



(16) 



where the dimensionless coefficient A has been intro- 
duced to retain the notation of the literature 30 . When 
k — ► oo, fxc(k) has the following asymptotic expansion 
(up to k~ 2 order) (see Appendix 



4^ 4tt 

rfjp A, 



(17) 



where C and B are also dimensionless quantities. 

The coefficients A, B, C naturally depend on the in- 
teraction chosen. A is given by the compressibility sum 
rule (see, e.g., Refs. Fl7l which is valid for any electron- 
electron interaction. C is related to the reduced second 
moment 82 of the momentum distribution. The expres- 
sion of B involves, among other things, the second and 
fourth moments 82 and 84. 

For the Coulomb interaction, ^4 cou i is given by 



.4, 



coul 



k F d 2 {ne xc 



-in 



dn 2 



(18) 



with the exchange-correlation energy per particle e X c.co\i\ 
taken from the usual VWN parametrization*!^. The coef- 
ficients C C oui and -Bcoui have been calculated by Holasi^S 
[see Appendix |0 Eq. (|C25f) ] and read 



a 



coul 



kpS2,cou\ 

5u> 2 



(19) 
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FIG. 11: Coefficient Seoul with the Coulomb interaction: GZ 
(solid curve) and MCS— (dotted curve) parametrizations, 
RPA (short-dashed curve) and Go Wo (long-dashed curve) cal- 
culations. 



where lu p = \/4irn is the plasma frequency, <? C oui(OHs the 
on-top pair-distribution function taken from Ref . 1 761 and 
82, coul is given by the virial formula of Eq. I|A5|I . 

The different coefficients -B C oui resulting from the GZ, 
RPA and GqWq parametrizations of e^coui [Eq. H-Bljl ] are 
compared in Fig. 1111 The parametrization of Moroni, 
Ceperley and Senatore (MCS)Si is also reported. The 
Go Wo calculation gives a reasonable estimation of S CO ui) 
in agreement with the GZ and MCS parametrizations. 

The coefficients A cou \, B cou i and C cou i have been cal- 
culated (using the G0W0 parametrization of #4, C oul) an d 
fitted to simple analytical formulas for convenience in 
Appendix [DJ 

For the long-range erf interaction, the coefficient A eI { 

is 



^erf 



k% 9 2 « C)Crf ) 
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dn 2 



(21) 



where ez. f is the exchange-correlation energy per par- 
ticle associated to the erf interaction taken from Ref. [27]. 
The calculation of Holas is generalized in Appendix IQfor 
the coefficients C er f and B cr f [Eq. (|C26() ]: they write 



k 4 8^ 
5uj 2 



(22) 



and 



B» r 

erf 



2/j, 



ut xc,crf 



3irn d/i 



ik F 



1 9 it 4 8^ 
LZKp 4 crf 

35^2 



(25L rf + (-5L rf ) 2 ) 



25lo 2 



(23) 



and 



Seoul (0)) 



I2kp8^ 



coul 



35^2 



4fc|,(2c> 2 , coul + (<5 2 .coul } ) 



25w 2 



(20) 



where 8% eI r is given by the virial formula of Eq. I|A6|I . The 
three coefficients have the correct Coulombic limit when 
H — > 00: A^°° = A couh B£j 

Ccoui ■ 

The coefficients A^ r{ , B^ f and C^ rf have been cal- 
culated (using the GqWq parametrization for 8% elf 



= B coul and C^°° = 
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FIG. 12: Static exchange-correlation kernel f X c,coui(k) with 
the Coulomb interaction for r s = 5: simple interpolation 
formula [Eq. 1251 . dashed curve], MCS^ (solid curve) and 
CDOP 31 (dotted curve) parametrizations. 



[Eq. (Bat ]) an d fitted to analytical formulas in Ap- 
pendix [DJ 

For both the Coulomb and long-range erf interactions, 
fxc(k) is then approximated over the whole range of k by 
the simple following stepwise interpolation between the 
two limiting behaviors at k — ► and k — > oo 



fxc(k) — 




fxc(0) >/£(*) 
/*c(0) < /£(*). 



(24) 



The junction point between the short and long wave vec- 
tor regions where f xc (0) — fSc(k) is located at k = y^fk F 
with 7 = B/(A — C), and Eq. 124|) can be re- written in 
the more compact form 



/-<*> - 
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B 



C- A 



'(fc- ^jk F ). 



(25) 



In the case of the Coulomb interaction and for r s = 5, 
Fig. E| compares interpolation of Eq. (I25|l with the 
parametrization of MCS [Eq. (7) of Rcf. 30 with ra = 8] 
and the parametrization of Corradini, Del Sole, Onida 
and Palummo (CDOP)2i based on the same QMC data. 
The three curves are in overall agreement, the main dis- 
crepancy arising around k ~ 2k F where the presence of a 
singularity or not is still a matter of debate22*22i2iiS. Of 
course, our simple model cannot be trusted for k « 2k F . 

Taking the inverse Fourier transform of Eq. I|25l) . we 
obtain the static exchange-correlation kernel for both the 
Coulomb and long-range erf interactions in real space 



4ttC 



fxc(r) = — p~ 5 ( r ) - — 



1 Si(y^/k F r) 



2(C — A) _ _ _ 

[sm(y/jk F r) - ^/jk F r cos(y^fk F r)] , (26) 



nr 3 kp 



where Si(x) — sm{t)/tdt is the sine integral function. 
In Eq. 1(26(1 . the second term gives the major contribution 



to fxcif) for r ^ 0, the third term being only a small 
correction. 

Fig. 1 131 shows the static exchange-correlation kernel in 
reciprocal and real spaces for r s — 5 with the erf inter- 
action for a series of interaction parameters fj,. f£ c orf is 
naturally reduced when the interaction decreases. More 
precisely, in reciprocal space, f^ cei {(k) is flatten and the 
junction point is slightly shifted toward larger values of k, 
while, in real space, the spatial extension of f£ c { (r) de- 
creases. The kernel f£ { thus becomes more and more 
local when fi decreases. Symmetrically, the complement 
kernel /£ c erf = / KC , cou i - /£ 
more local when fj, increases (not shown) 



M c orf also becomes more and 



In the limit of the Coulomb interaction, the kernel 
of Eq. I|2t)|l displays Friedel-like long-range oscillations 
in real space of quasiperiodicity -Jjk F » 2k F coming 
mainly from the second term in Eq. I|26|) . These os- 
cillations connected with the behavior of the kernel in 
reciprocal space around k w 2k F constitute an usual fea- 
ture of models for the exchange-correlation kernel of the 
uniform electron gas. In our model, these oscillations 
are likely to be exaggerated due to the sharp cut-off in 
Eq. (|25|l and their relevance for inhomogeneous systems is 
questionable. However, the oscillations in our model for 
fxc erf ( r ) ra pidly spread out and vanish when fi decreases, 
i.e. when the interaction is reduced, (see Fig. 113(1 so that 
this problem disappears. Symmetrically, the oscillations 
in the the complement kernel f£ 01 .((r) also disappears 
when (x increases. This results are consistent with the 
common intuition that the decrease of the interaction re- 
duces the structure of the exchange-correlation kernel. 
Note that in the uniform electron gas the interaction can 
also be reduced by decreasing r s so that an equivalent 
way to look at fl^ cer {{r) in the limit of a small interac- 
tion, (i — ► for a fixed r s , is to look at the high-density 
limit, r s — * 0, for a fixed 



IV. CORRELATION ENERGY FROM STATIC 
EXCHANGE-CORRELATION KERNEL 



As an example of the usefulness of the static exchange- 
correlation kernel, we follow Lein, Gross and Perdew21 
and compute the correlation energy of the uniform elec- 
tron gas from it. 

Indeed, the correlation energy per particle can be ex- 
actly deduced in principle from the dynamical exchange- 
correlation kernel using the ACFD approach. For com- 
parison purposes, we begin by the standard case of the 
Coulomb interaction for which the expression is (see, e.g., 
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rs = 5 . 



f£c,erf (*) 
(a.u.) 




rs = 5 . 





-0 . 02 

f xc,erf( r ) -0.04 

(a.u.) 

-0 . 06 
-0 . 08 




FIG. 13: Static exchange-correlation kernel f£ c orf in reciprocal [Eq. |25p ] and direct [Eq. I|26|l ] spaces for r s — 5 with the erf 
interaction (dashed curves) for fj, = 0.2, 0.3, 0.5 and 1.0. For comparison, the case of the Coulomb interaction (solid curve) is 
also shown. 



Ref.Ei 



_ACFD 



— 1 

27m Jq 



dto- 



Vcoui(q) 





U coul(?) + /xc,coulfo iw ) 




1 




U coul(9)+/xc,coul(9^ W ) 





(27) 



with the scaled interaction v£ oul (q) = Xv cou \(q) and the 
associated exchange-correlation kernel f£ c coul (<Z, u>) re- 
lated to /xc,coul(<7)k') by a simple scaling relational 

/icou! N (9, w) = A-V.ccouifn/A 3 ] (g/A, ^/A 2 ). (28) 

We can thus use Eq. I|27|l to estimate the relevance of our 
simple model of the static exchange-correlation kernel to 
obtain correlation energies, even though in theory the 
frequency-dependence is also needed to obtain the exact 
correlation energy. Fig. 1141 reports the error in the corre- 
lation obtained with Eq. (j2ZJ>, Ae C:COU i = ~ £ c,coui, 
using several approximations for / xc ,coul- RPA (/ xc ,coul = 
0), MCS and CDOP parametrizations and the interpo- 
lation formula of Eq. 1251) . The expansion of the exact 
correlation energy for r s — > (see, e.g., Refs. l7§l79ll80l) is 
£ c ,coui = C In rv, + Ci + 0r s \nr s with C Q = (1 - In2)/7T 2 
and d = -0.0469205 while in the RPA it is e^i = 
C In r s + Cf PA + 0r s lar a with C? PA = -0.071100. The 
error in the RPA correlation energy at r s — > is due to 
the difference Cf = C x - Cf PA = 0.024179 correspond- 
ing to the contribution of exchange diagrams. When a 
static exchange-correlation kernel is used in Eq. (|27|) , the 
correlation energy is greatly improved compared to the 
RPA. The exchange-correlation kernel of Eq. 1|25|1 gives 
results very close to those given by the MCS and CDOP 
parametrizations, validating our simple interpolation. 

We now generalize the ACFD approach to the case of 
the long-range erf interaction. In this case, the long-range 




. 005 



-0.005 

AEc.coul 

(a.u.) -0.01 
-0.015 
-0.02 
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r s (a.u.) 

FIG. 14: Error in the correlation energy Ae c>cou i = ejr^ju? — 
£c,coul for the Coulomb interaction [Eq. 1271 1 with differ- 
ent exchange-correlation kernels: RPA [f xc = 0, short- 
dashed curve], MCS^ (solid curve) and CDOP^i (dotted 
curve) parametrizations and the simple interpolation formula 
[Eq. 1251 . long-dashed curve]. 



correlation energy per particle is given by 

<iq 



_/i,ACFD _ 



-1 

27T71 



(27T) 



dto- 





U erf(9)+/xcJrf(2>H 




1 




«rf (?) + /x<£rf(?. iuJ ) 





(29) 

with the scaled interaction v^(q) = \v^ lf (q) and the 
associated exchange-correlation kernel f^ eeti (<?> uj). The 
scaling relation of Eq. (|28ll is easily generalized to the erf 
interaction (see also Ref. l8lh 



A-VxterfKA^Cg/A^/A 2 ). (30) 



The error in the long-range correlation energy obtained 



_jU,ACPD 
"c,erf 



<erf> in the RPA 



with Eq. ((23, Ae£ erf 
(fxc erf = 0) an d with the exchange-correlation kernel of 
Eq. (|25f) are represented in Fig. ^| for \i — 1. The reduc- 
tion of the interaction mainly acts at small r s decreasing 
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^=1.000 
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A< erf 
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FIG. 15: Error in the correlation energy Ae^ erf = £^'^.f FT> — 
£ c erf f° r the erf interaction [Eq. I29H 1 with different exchange- 
correlation kernels: RPA [f xc = 0, short-dashed curve], the 
simple interpolation formula [Eq. 1251 . long-dashed curve], for 
fi = l. 



APPENDIX A: KINETIC ENERGY FROM 
VIRIAL THEOREM 

In this appendix, using the virial theorem, we give the 
expression of 62 [Eq. for the Coulomb and erf inter- 
actions. 



the error given by the RPA. As for the Coulomb case, the 
exchange-correlation kernel of Eq. (|25l) enables to correct 
the RPA for larger values of r s . 



V. SUMMARY AND CONCLUSION 

We have given a simple approximate expression in real 
and reciprocal spaces for static exchange-correlation ker- 
nel f£ orf of a uniform electron gas interacting with the 
long-range part only of the Coulomb interaction [Eq. i|25|) 
and Eq. (|26j . This expression interpolates between the 
exact asymptotic behaviors of the kernel at small and 
large wave vectors which in turn requires, among other 
thing, information from the momentum distribution of 
the uniform electron gas with the same interaction that 
we have calculated in the GqWq approximation. In the 
limit of the Coulomb interaction (p, — > 00), the proposed 
exchange-correlation kernel /^T^Jf = /i C , C oui is closed to 
other parametrizations of the literature based on QMC 
data. As a matter of course, the complement kernel 
fxc,er{ = fxc : cou\ - /^ C;0rf associated to the short-range 
part of the Coulomb interaction is also obtained in this 
work. 

In the same way that the static exchange-correlation 
kernel of the Coulombic uniform electron gas is an essen- 
tial ingredient in ACFD and DFT approaches of the elec- 
tronic correlation problem, the modified kernels f£ f 



' c orf can be used in alternative ACFD and DFT 



and /£ 

approaches based on a separate treatment of long-range 
and short-range correlations. As a simple illustration, 
the long-range correlation energy of the uniform electron 
gas has been calculated from f£ c erf . In comparison to the 
Coulombic kernel f xc , cou i, the modified kernels f£ e orf and 
fxc orf are more local which may improve their transfer- 
ability to inhomogeneous systems and therefore facilitate 
their use in applications. We are currently investigating 
such an issue. 



1. The virial theorem 

The virial theorem for the uniform electron gas with 
the Coulomb interaction writes 82,83 



2^COUl ) + v cou i{r s ) 



de cou i{r s ) 
dr. 



(Al) 



where iconics) = to(r s ) + t CtCOU i(r s ) is the interacting 
kinetic energy per particle which can be decomposed 
into a non- interacting contribution to(r s ) and a corre- 
lation contribution t CjCOU i(r s ), v cou \(r s ) is the electron- 
electron repulsion energy per particle and e cou \(r s ) = 
iconics) + Wcoui^s) = to{r s ) + s XC} coui{r s ) is the total 
energy per particle which can be decomposed into the 
non-interacting kinetic energy to(r s ) and the exchange- 
correlation energy per particle s xc ^ cou i(r s ). Eliminating 
u cou i(r s ) for £ xc ,cou\( r s) in Eq- <|A1|) and using the virial 
theorem for the non-interacting electron gas, 2to(r s ) = 
—r s dto(r s )/dr s , and the virial relation for exchange only, 
£x,coul(^) = —r s ds XtCOU i(r s )/dr s , we arrive at the virial 
theorem for the correlation part (see also Ref. IHEh 

tc,coul{ r s) F £c,coul(?*s) = r s ; ■ (A2) 

dr s 

For the erf interaction, the virial theorem of Eq. (|A1|) 
is generalized as (see Ref. l8lj) 

«+< rf W-^ = -^. (A3) 



Following the same steps as for the Coulomb interaction, 
it can be shown that the virial theorem for the correlation 
contributions reads now 



*c.erf( r s)+ £ c.erf( r s)-M 



9 <orf( r s) _ 9e c,erf( r s) 



Of 1 



dr s 



(A4) 
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2. Expression of 82 

Using t c , cou i(r s ) = -d(r s e c ^ cou \(r s )) / dr s coming from 
Eq. (|A2|) and t (r s ) — 3k F /10, the expression of 5 2 . C oui 
for the Coulomb interaction writes 



t. 



c.coul 



- , 2,coul 



to(r s ) 

10 d(r s e C)COU i(r s )) 



3k p 



(A5) 



For the erf interaction, Eq. (|A4|) gives t^ eT{ (r s ) = 
-d(r s e% exf (r s )) /dr s + fide% eli (r s ) /dfi and the expression 



0f 5 2,crf is 



°2,erf — 



*o(r.) 

JL0_ /_ 9(r s < erf (r s )) 
3fc| 1 



<3r s 



t ^ £ c,erf( r ») 

dfj, 



(A6) 



APPENDIX B: ANALYTICAL 
PARAMETRIZATION OF <5 4 

In this appendix, we give analytical paramctrizations 
for the reduced moment 64 of the momentum distribution 
[Eq. Q14|)] calculated by several methods for the Coulomb 
and erf interactions. 



Coulomb interaction 



For the Coulomb interaction, 8^ cou \ is parametrized as 



&31 


= -0.690727, e 32 


= 0.393525, 


e 33 = 


3.398631, 




= 0.231529, e 42 


= 0.885698, 


e43 = 


5.999882, 




= 0.001233, e 52 = 


= -0.061940, 


e53 = 


2.887302, 


eei 


= 0.001083, e 62 = 


= -0.000060, 


e63 = 


0.391712. 


Those for the GqWq calculation are 


E31 = - 


-0.525472, 


e32 


= 0.417720, e 33 


= 5.102281, 


e 4i = 


0.223039, 


e42 


= 0.316492, e 43 


= 5.984490, 


e-51 


0.002486, 


e52 


= -0.000770, e 53 


= 0.265086, 


1 IV = - 


-0.004998, 


e62 


= -0.003175, e 63 = 


1.191797. 







APPENDIX C: ASYMPTOTIC EXPANSION OF 
THE STATIC EXCHANGE-CORRELATION 
KERNEL 

In this appendix, following the procedure proposed by 
HolasZ 4 * 7 ^, we derive the asymptotic expansion of the 
static exchange-correlation kernel f X c(k) of the uniform 
electron gas for k — > 00 for the erf interaction. For 
comparison, results for the Coulomb interaction are also 
given. Contrary to the rest of this work and except for 
Eqs. I|C25|) and (|C26|) . momentum are expressed in units 
of kp , energies and frequencies in units of k F throughout 
this appendix. 



Frequency moments of the linear response 
function 



Knowledge of the frequency moments of the frequency- 
dependent linear response function x(k,u>) is useful to 
study the asymptotic behavior of the static response 
function for k — > 00. These moments are defined as 



Mi{k) 



Imx(fc, w)w dw. 



(CI) 



^4,coul 



2 = 3 



(Bl) 



where x = y^rv A least-square fit of the result of the GZ 
parametrization gives d 3 = 0.271191, d A = -0.009998, 
d 5 = -0.036383, d 6 = 0.006706. For the RPA cal- 
culation we obtain d 3 = 0.093623, d 4 = 0.194288, 
d 5 = 0.051445, d 6 = 0.005449, and for the G a W cal- 
culation, d 3 = 0.126362, d 4 = 0.001428, d 5 = 0.014278, 
d 6 = -0.004522. 



erf interaction 



The zeroth, first and third moments are well- 
known 74 ! 75 . 84 



M 3 (k) 



M {k) = -2nS(k), 



Mi(fc) = -nfc 2 , 



- nk 1 { — + 2fc 2 (E K ) + nk 2 v{k) 



-47rnG FV (/c) , 



(C2) 



(C3) 



(C4) 



For the erf interaction, we take the parametrization 
< erf = X>^ (B2) 



i=3 



with e% = (enfi + e t2 [i 2 )/{l + e i3 ^ 2 ). The opti- 
mal parameters reproducing the RPA calculation are 



where S(k) is the static structure factor, Ek — 
(1/N) X)kcr e (^) c k(T C kcr is the kinetic energy per parti- 
cle operator, (• • • ) means expectation value between the 
interacting ground state, and G PV (/c) is the "Pathak- 
Vashishta"^ (or "Niklasson"^ 5 .) local-field factor 



G^(k)=G^(k) + Gr(k) 



(C5) 
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with 



and 



If dp (k • p) 2 

(2tt) 3 ~ P 

dp (k-p) 2 v(p) 



(S(p) - 1) , (C6) 



(2tt) 3 k 2 



4tt 



(5(|k- 



The asymptotic expansion of M (fc) as fc 
by that of 

M (k) = -2n(l + £ + (i 



■Pl)-1). 

(C7) 
oo is given 



(C8) 



For the Coulomb interaction, C is directly related to the 
cusp condition: C CO ui = — 2uj 2 g'(0) — —2uj 2 g(0) where 

uj p = \J Aim is the plasma frequency, 17(0) and g'(0) are 
the on-top values of the pair-distribution function and 
of its spherical-average derivative. For the cusp-less erf 
interaction, C vanishes: C er f = 0. 

The asymptotic expansion of M^{k) writes 

M 3 (fc) = -^(l + ^ { E K ) + ±nv(k) 



k 



(C9) 



where the limit G PV (fc — ► 00) is determined as follows. 
G„(fc) is actually independent of k and writes 



11/" dp p 2 v(p) 



(S(p)-1). 



in J (2tt) 3 4tt 
which gives, for the Coulomb interaction, 



(CIO) 



G 



PV 



1 1 

3 n 
1 
3 



dp 



(2tt) 3 

(ffcoul(O) - 1) 



(Seoul (p) ^ 1) 



(Cll) 



and, for the erf interaction, 



G 



1 1 

3n 

67m 



M 3 dS xc.m-l 

37m 9/i 



dp 

(2^r)3 

dp <9< rf (p) 
(2tt) 3 dp 



c -pW) ^(p) _ i) 



(C12) 



where the Hellmann-Feynman theorem has been used for 
the last line, and e^ c erf is the exchange-correlation en- 
ergy per particle with the erf interaction. 

After a trivial variable transformation, the limit of 
00 is easily seen to be 



G ^ as k 



✓-.PV 



(fc — ► Oo) = 



1 f dp (k-p) 2 «(|k- P |) 



n J (2tt) 3 
x (S{p) - 1) , 



47T 



(C13) 



which gives, for the Coulomb interaction, 

GbJ ou i(fc -> 00) = ~- [ 7^3 (5coul(p) - 1) 



(ffcoul(O) - 1) , 



(C14) 



PV,M/i. _ _I / rf P -(k-p) 2 /(4 M 2 ) 



and, for the erf interaction 

- - n</ (27r) 3- 

x(s£f(p)-i) 

= -e-W)(^(0)-l), (C15) 

where <?e rf (0) is the on-top pair-distribution function as- 
sociated to the erf interaction. 



2. Static linear response function 

The static linear response function x(fc,0) can be ex- 
pressed by the spectral representation 



X(k,0) 



lmx{k,ui) 



du>, 



(C16) 



i.e., we simply have x(fc,0) = M_i(fc). Using this re- 
sult, HolasZ^iS showed that the asymptotic expansions 
for large k of M (k), M x {k) and M 3 (k) [Eqs. JUS)), (IH3l) 
and (|C9|) ] are sufficient to determine the asymptotic ex- 
pansion of x(k,0) up to fc~ 4 order 



x(M) 



—An 



8,- ,1 / 1 
1+ 3™F + ° IF 



(C17) 



He then determined the following terms in the asymptotic 
expansion of x(fc, 0) by (not rigorously) inferring them as 
follows. The following term in k~ 3 in the square-bracket 
of Eq. (|C17() is set to zero by a simple argument involv- 
ing the continuity and the likely non-oscillating behavior 
of the moments Mi(k) with respect to I. The next term 
in fc~ 4 is inferred from the asymptotic expansion of the 
static response function calculated to first-order with re- 
spect to the electron-electron interaction. We follow the 
same procedure. 

The (zeroth-order) Lindhard static response function 
Xo(k, 0) is well-known. Its asymptotic expansion for large 
k writes 



Xo(fc,0) 



—An 



Y {Ek)o ¥ 



(cis) 



where E 2 K = (1/JV) J2ka £ ( k ) 2c i a c ^ and ("Oo means 
expectation value between the non-interacting ground 
state. 

The asymptotic expansion of the first-order correction 
to the static response function xi(k,0) has been calcu- 
lated for the Coulomb interaction by Geldart and Tay- 
loi-^Ml (see also Ref. Ia3j) . This result is conveniently 
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generalized to an arbitrary interaction as 



xi(M) 



—An 
IF 



-AojI f 2G PV 



fc 2 u(fc) 



HF 



4tt 



with 



G PV 



HF 



1 1 

3 n 



fc 4 



dp p 2 v(p) 



1 

F 

(C19) 



(5hf(p) ~ 1) , (C20) 



(2tt) 3 4tt 

and the Hartree-Fock (HF) static structure factor 
2 f dk 



Shf(p) = 1 - 



(2tt) 3 



no(A;)no(|k + p|). (C21) 



For the 



^a.HF.coul 



Coulomb interaction, we simply have 
= (.9hf,coui(0) - l)/3 where 

<?HF,coul 

(0) = 1/2 

is the HF on-top pair-distribution function. For the erf 
interaction, G^ >erf = ^ / \3im)de^ cii /dv where e» f 
is the exchange energy per particle associated to this 
modified interaction. 

Following HolasZ 4 *^, we then infer the following term 
in fc~ 4 in the square-bracket of Eq. (|C17|1 from the corre- 
sponding term in the expansion of Xo(k, 0)+xi(fc, 0) with 
the substitutions (E 2 K ) -> (E 2 K ) and G P ^ F -> G PV . We 
thus finally obtain 



The term explicitly depending on v(k) in the asymp- 
totic expansion of 0) [Eq. l|C22|l ] exactly cancels 
with the Hartree kernel v(k) in the expression of fxcik) 
[Eq. 1C23I) ]. Therefore, the asymptotic expansion of 
fxcik) [Eq. i|C24|) ] do has the same form for the Coulomb 
and erf interactions. Note however that the local field 
factor G(k) = —v{k)f xc (k) involves explicitly the inter- 
action and has therefore different asymptotic behaviors 
for the Coulomb and erf interactions. 

Coming back to atomic units now and specializing to 
the Coulomb interaction, Eq. 1C24() becomes 



IEC, COUl 



(k) 



-47rfc|(5 2 ,caul 



35w! 



44 



4tt (2, 

F U (1 



Seoul (0)) 



(2(52, coul + (<^2,coul) 2 ) 

"25^ 



1 

F 

(C25) 



while for the erf interaction, we have 



/xcerf(^) 



5^2 



An 

F 



2/Z 3 < 6 

37rn 9/i 



35^2 



44(2^, err 



(<*U) 2 : 



25^2 



(C26) 



X (fc,0) = 



-4n 



1 



-4w; 2G 



-.PV 



k 2 v(k) 

47T 



1 

F 



(C22) 



Observe that the form of the asymptotic expansion of 
x(k, 0) [Eq. (|C22Jl ] does depend explicitly on the interac- 
tion v(k) and is therefore different for the Coulomb and 
erf interactions. 



APPENDIX D: ANALYTICAL 
PARAMETRIZATIONS OF A, B AND C 

To facilitate the evaluation of the static exchange- 
correlation kernel for the Coulomb and erf interaction, 
we give analytical parametrizations of the coefficients A, 
B and G determining its limiting behaviors for small and 
large k (see Sec. 



3. Static exchange-correlation kernel 

The static exchange-correlation kernel f X c(k) writes 
[see Eq. ifTKJl] 



fxc(k) = Xo(k, 0)- 1 - X (k, 0)- 1 - v(k), 



(C23) 



1. Coulomb interaction 

For the Coulomb interaction, the coefficient A 
[Eq. l(T5)l ] is parametrized as follows 



A-coul — / &iX 



(Dl) 



and its asymptotic expansions for large k is therefore 
determined from the expansions of Xo(fc>0) and %(fc,0) 
[Eqs. (|C18(I and (|C22|) ]. Introducing the quantities 5 2 = 
{{E K ) - (E K )o)/(E K ) and 5, = ((E 2 K ) - (E 2 K ) )/(E 2 K ) 
with (E K ) = 3/10 and (E%) = 3/28, we obtain (with 
k in units of fcp and uj p in units of k 2 F ) 



fxcik) 



-Att5 2 

" p 



f 1 



12S 4 4(2<5 2 + 51) 



35u;2 



25^2 



(C24) 



with x = y/r^ and the fitted parameters ao = 0.250019, 
ai = -0.000162, a 2 = 0.013441, a 3 = -0.003591, a 4 = 
0.000380, a 5 = 0.000002, a 6 = -0.000003. 

For the coefficient B (Eq. I20fl , we adopt the analytical 
form of Moroni, Ceperley and Senatore 3 ^ incorporating 
the correct limit -B(r s — > 0) = 1/3 



-^coul — 



1 + b\x + b 2 x 3 
3 + 63X + 64a; 3 



(D2) 



but the parameters are refitted with the parametrization 
of 8 A [Eq. ljBl|l ] with the G VF data: 61 = 0.721543, 
62 = 0.317320, 63 = -0.133379, 64 = 0.269494. 
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The coefficient C [Eq. (|19fl ] is parametrized as 

Ccoui = 1 „ — T j (D3) 



with ci = 0.002127, c 2 = 0.169597, c 3 = 0.450771, 
c 4 = -0.023265, c 5 = 0.001855, c 6 = -0.000069, /i = 
7.062604, / 2 = 8.589773, / 3 = 2.747407, / 4 = 0.648920. 



2. erf interaction 

For the erf interaction, A^ r{ [Eq. (|2*T|l ] is represented 
by the analytical parametrization 



v^6 i 

crf ~ l + /i 4 x 4 



(D4) 



with gi — f/ii/i + gi2^ 2 and /14 = 341/i + /142/i 2 . The fit- 
ted parameters are 321 = -0.029315, 522 = -0.000927, 
.931 = 0.061867, g 3 2 = 0.0 1 09 70, 341 = -0.053761, 
342 = 0.078580, g 5 i = 0.012970, 352 = 0.014669, 
36i = -0.001232, 3 6 2 = -0.000891, hu = -0.025963, 
h 42 = 0.389673. 

The coefficient B^ r{ [Eq. J33J)] is parametrized as 



Bt 



E3 - i 
i=Q 3i x 

1 + k\x + k 3 x 3 ' 



(D5) 



where ]\ = (jafi + j l2 fi 2 )/(l + 3i3^ 2 ) and h = [kafi + 
ki2H 2 ) I (l+kisfi 2 ) . Using the parametrization of Eq. (|B2|) 
for S% in the Go Wo approximation, the fitted param- 
eters obtained are j i = 0.010533, j 2 = -0.002640, 
= -0.143455, ji2 = 0.046302, 
= -0.415043, j 2 2 = 0.194149, 
= 0.164085, J32 = 0.925083, 
= 0.291633, k 12 = 0.102905, 
0.014324, k 32 = 0.7 1 44 49, h 



303 

jn 

323 
333 

ki 3 



= 0.314403, jn 
: 0.014315, j 21 
-- 0.078849, j 31 
= 0.491612, k n 
0.013551, k 31 = 



33 



0.473130. 



Finally, the coefficient C^. f [Eq. I|22|l ] is parametrized 



as 



^erf 



1 + m^x 4 



(D6) 



' 2 )/(l 



where li 
(1TI41H + m 42 ^ 
are l n = -0.039662, 
/ 2i = 0.187782, Z 22 
I31 = -0.270823, Z 32 
l 41 = 0.133824, l 42 = 



^i2ju 2 )/ (1 + li3fJ- 2 ) and m 4 = 
w 4 3/i 2 ). The fitted parameters 
l u = 0.002346, l 13 = 0.405492, 
= 0.066673, l 23 = 0.759503, 
= -0.083482, Z 33 = 0.511806, 
0.061658, Z 43 = 0.348085, /51 = 



-0.027392, l 52 = 0.02 4 3 84, l 53 = 0.26 1 739, l 61 
0.001877, l e2 = -0.005601, l 63 = 0.35 1 2 72, m 4 i 
-0.090916, m 42 = 1.238974, m 43 = 0.25 1 286. 
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